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Abstract 

We consider a cluster growth model on called internal diffusion limited aggre- 
gation (internal DLA). In this model, random walks start at the origin, one at a time, 
and stop moving when reaching a site not occupied by previous walks. It is known 
that the asymptotic shape of the cluster is spherical. When dimension is 2 or more, we 
prove that fluctuations with respect to a sphere are at most a power of the logarithm 
of its radius in dimension d > 2. In so doing, we introduce a closely related cluster 
growth model, that we call the flashing process, whose fluctuations are controlled easily 
and accurately. This process is coupled to internal DLA to yield the desired bound. 
Part of our proof adapts the approach of Lawler, Bramson and Griffeath, on another 
space scale, and uses a sharp estimate (written by Blachere in our Appendix) on the 
expected time spent by a random walk inside an annulus. 
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1 Introduction 



The internal DLA cluster of volume N, say A(N), is obtained inductively as follows. Initially, 
we assume that the explored region is empty, that is A(0) = 0. Then, consider N independent 
discrete-time random walks Si, ... , Sn starting from 0. Assume A(k — 1) is obtained, and 
define 

r fe = inf{t>0: S k {t) £ A(k - 1)} , and A(k) = A(k - 1) U {S k (r k )}. 
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In such a particle system, we call explorers the particles. We say that the k-th explorer is 
settled on S^t^) after time r^, and is unsettled before time t&. The cluster A(N) consists of 
the positions of the N settled explorers. 

The mathematical model of internal DLA was introduced first in the chemical physics 
literature by Meakin and Deutch [13]. There are many industrial processes that look like 
internal DLA (see the nice review paper [7]). The most important seems to be electropol- 
ishing, defined as the improvement of surface finish of a metal effected my making it anodic 
in an appropriate solution. There are actually two distinct industrial processes (i) anodic 
levelling or smoothing which corresponds to the elimination of surface roughness of height 
larger than 1 micron, and (ii) anodic brightening which refers to elimination of surface defects 
which are protruding by less than 1 micron. The latter phenomenon requires an understand- 
ing of atom removal from a crystal lattice. It was noted in [T3] that, at a qualitative level, 
the model produces smooth clusters, and the authors wrote "it is also of some fundamental 
significance to know just how smooth a surface formed by diffusion limited processes may 



Diaconis and Fulton [2] introduced internal DLA in mathematics. They allowed explorers 
to start on distinct sites, and showed that the law of the cluster was invariant under per- 
mutation of the order in which explorers were launched. This invariance, named the abelian 
property, was central in their motivation. They treat among other things, the special one 
dimensional case. 

In dimension two or more, Lawler, Bramson and Griffeath [10J prove that in order to 
cover, without holes, a sphere of radius n, we need about the number of sites of Z d contained 
in this sphere. In other words, the asymptotic shape of the cluster is a sphere. Then, Lawler 
in [9] shows subdiffusive fluctuations. The latter result is formulated in terms of inner and 
outer errors, which we now introduce with some notation. We denote with || • || the euclidean 
norm on M d . For any x in IR d and r in M, set 



For A C Z d , |A| denotes the number of sites in A. The inner error Si(n) is such that 



be". 



B(x,r) = {y E R d : \\y - x\\ < r} and M(x,r) = B(x,r) (1 Z d . 



n-S^n) =sup{r > : 1(0, r) C A(|1(0, n)\)} . 



Also, the outer error 5o{n) is such that 



n + 5 (n) = inf {r > : A(|B(0,n)|) C 1(0, r)} . 



The main result of [9] reads as follows. 



Theorem 1.1 [Lawler] Assume d > 2. Then, 



P (3n(w) 



Vn > n{u) Sj(n) < n 



/3 log(n) 2 ) 



1- 



(1.1) 



and 



P (3n(w) 



Mn > n{oS) S (n) < n 1/3 log(n) 4 ) 



1. 



(1.2) 
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Since Lawler's paper, published 15 years ago, no improvement of these estimates was achi- 
eved, but it is believed that fluctuations are on a much smaller scale than n 1//3 . Moreover, 



(1.1) are almost sure upper bounds on errors, and no lower bound on the inner or outer 
error has been established. Computer simulations [HJ E] suggest indeed that fluctuations 
are logarithmic. In addition, Levine and Peres studied a deterministic analogue of internal 
DLA, the rotor-router model, introduced by J.Propp [6j. They bound, in [12], the inner 
error <5/(n) by log(n), and the outer error So(n) by n 1 ^ 1 ^. 

We now introduce a family of cluster growth models for which a control of the fluctuations 
of the cluster shape is easily obtained. These growth models are built so that the asymptotic 
shape is spherical, but still they exhibit a large diversity of fluctuations parametrized by a 
certain width ranging from a large constant to a power 1/3 of the radius of the asymptotic 
sphere. Moreover, all these clusters are coupled to internal DLA, and, as a consequence, we 
obtain logarithmic bounds on the fluctuations for internal DLA. We generalize internal DLA 
by allowing explorers to settle only at some special times. Thus, each explorer i is associated 
with a collection of times {o^, k G N} and 

r* = infK fe : S^) # A*{i - 1)} , and A*(i)=^(»-l)|J{» 

The internal DLA is recovered as we choose cr^ = k for all i — 1, . . . , N and k EN. We call 
{o^fc, k G N} the flashing times associated to the i-th explorer, and {^(er^), k G N} its 
flashing positions. 

In this paper, we consider stopping times of a special form, linked with the spherical 
nature of the internal DLA cluster. An illustration with one flashing explorer's trajectory is 
made in Figure [T] 




Figure 1: Cell decomposition, and flashing positions as stars. 

The precise definition of the flashing times requires additional notation, which we post- 
pone to Section [3j We describe here key features of flashing processes. We first choose a 
sequence of widths, say H = {h n , n G N}, and then partition Z d into concentric shells 
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{S n , n e N}, whose respective widths are {2h n , n e N}. Each shell is in turn partitioned 
into cells, which are brick-like domain, of side length equal to the width of the shell. The 
flashing times are chosen such that (i) an explorer flashes at most once in each shell, (ii) the 
flashing position, in a shell, is essentially uniform over the cell an explorer first hits upon 
entering the shell, and (iii) when an explorer leaves a shell, it cannot afterward flash in it. 

For a given sequence H, we call the process just described the ^-flashing process. Note 
that feature (ii) is the seed of a deep difference with internal DLA. The mechanism of covering 
a cell, for the flashing process, is very much the same as completing an album in the classical 
coupon- collector process. Thus, we need of the order of Vlog(V) explorers to cover a cell 
of volume V. For internal DLA, with explorers started at the origin, we need only of order 
V explorers to cover a sphere of volume V as shown in [10J, and we believe that we need 
a number of explorers of order \C\ to cover a cell C, even if they start on the boundary of 
the cell. In addition, feature (ii) allows to localize the covering mechanism, in the sense that 
a particle entering a shell cannot flash outside the cell by which it entered into that shell. 
Finally, feature (iii) is essential for having a useful coupling between flashing and internal 
DLA processes. 

Lemma 1.2 Assume that N is an integer, and % is a sequence of positive integers. There 
is a coupling between the two processes, using the same trajectories Si, ... , Sn such that 

N N 

A(N) = {J{Si(T(i))}, and A*(N) = \J{Si(T*(i))}, (1.3) 

i=i i=i 

and T*(i) > T(i) for all i = 1, . . . , N. 



As a corollary of Lemma 1.2, we have the following useful result. 
Corollary 1.3 Under the hypotheses of the previous lemma, for k > 1, 

• if A*(N) C U j<k Sj, then A(N) C U j<k Sj, 

• tfU j<k Sj C A*(N), then U j<k Sj C A(N). 

An "H-flashing process, with hj > h for j > 0, and h a large constant, produces a cluster 
A*(N), for which we bound easily the inner error, 5j(n). Then, to bound the outer error, 
8q (n), we follow the approach of [9], though with a slightly simpler proof. 

Proposition 1.4 Assume that for j > 1, hj < hj + % < (1 + ^)hj, with a large Hq. For a 
positive constant A* d , we have 

P (3n(u) : > n(u) S*j(n) < A* d h(n) log(n)) = 1, (1.4) 

and, 

P (3n{u) : Vn > n(co) 5* (n) < A* d h{n) log 2 (n)) = 1, (1.5) 
where h{n) = max{h k G R : r k < n}. 
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Finally, we establish lower bound on the inner and outer error. 
Proposition 1.5 Assume that ho is large enough. Then, there is a constant a* d such that 
P (3n(w) : Vn > n{u) 5}(n) > a* d h(n) log(/i(n))) = 1, (1.6) 

and, 



P(3n(uj) : Vn > n(w) 5* (n) > a* d h(n)\og(n)) = 1. 



;i-7) 



Corollary 1.3 and Proposition 1.4, with the choice /i,- = /iq for all j > 0, imply the 



following result for internal DLA. 

Theorem 1.6 Assume d > 2. There is a positive constant Ad such that, 

P(3n(w) : \/n > n(u) <J/(n) < A d log(n)) = 1, 

and, 

P (3n(w) : Vn > < A d \og 2 (n)) = 1, 



;i.9) 



Our approach builds on the work of Lawler, Bramson and Griffeath [10J, which we first 
review. We then sketch our main ideas leading to logarithmic fluctuations for internal DLA. 



On previous bounds for internal DLA. We describe the approach of [TO], for estab- 
lishing the upper bound for the inner error. It is convenient to consider explorers starting 
outside the origin with initial configuration denoted rj. We denote also by A(A, rf) the cluster 
obtained from explorers initially on 77, with an explored region A C Z d . 

Now, for a site z G Z d , we call W(r), z) (resp. M(t], z)) the number of explorers (resp. of 
random walks) which visit z before settling. For an integer n, and 77 consisting of |B(0,n)| 
explorers at the origin, the authors of [TU] first write 

{1(0, r) <£A(Q,r))}c [j {W( V ,z)=0}. 

2SB(0,r) 

Then, they look for the largest value of r n (in terms of n) which guarantees that |B(0,r n )| x 
su PzeB(o,r„) P(W / ( ? 7) z ) = 0) De the term of a convergent series. 

The approach of [10] is based on the following observations, (i) If explorers would not 
settle, they would just be independent random walks; (ii) exactly one explorer occupies each 
site of the cluster. Thus, the following equality holds in law. 

W(rj, z) + M(A(0, r)),z)> M(r), z). 

Now, an observation of Diaconis and Fulton [2] is that we can realize the cluster by sending 
many exploration waves. Let us illustrate this observation with two waves. We first stop 
the explorers on the external boundary of a ball of radius R, say dM(0,R). We denote 
with Ar(0, rj) the cluster consisting of the positions of settled explorers, so that Ar(0, 77) C 
B(0, R). We denote by Cfl( r /) the configuration with stopped explorers on <9B(0,P). Then, 
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the second wave consists in launching the explorers of (nit]), with explored region Ar($,T]). 
In other words, we have an equality in law 

A($, V )=A R ($, V )UA(A R ($,r ] ),CR(v))- 
Moreover, if the index R refers only to explorers (or walks) of the first wave, then for 

z g 1(0, R) 

W R (r], z) + M R {A R {$, r)),z)> M R {rj, z). (1.10) 

The authors of [10J consider R = n and z G B(0,n). Since W(n, z) > W n (n,z), we have 
using (1.10), for any a > 

P (W{r), z) = 0)<P (M n {r), z)<a) + P (M n (B(0, n),z)>a). (1.11) 

We then look for sites z such that E[M n (r),z)] > a > E[M n (M(0,n), z)) (and rj = NS ). 
Note that M n (n, z) and M n (B(0, n), z) are sums of independent Bernoulli variables with well 
known large deviation estimates. If we set 

2a = E[M n (r],z)]+E[M n (M(0,n),z)], and fl n (z) = E[M n ( V , z)} - E[M n (M(0, n), z)], 

then, 

P(M n ( V ,z) <a)< exp (- {E[ ^lf \f \ < exp (- \ . (1.12) 

1 W ' V 2E[M n {r),z)] J' y \ 8E[M n {r),z))J 1 ; 

Lawler in [9] establishes that for z G B(0, n) 

E[M n (r], z)\ ~ n(n - \\z\\), and jx n (z) ~ (n - ||z||) 2 . 



Replacing these values in (1.12), the bound n— 1|2;|| > n 1 ^ 3 \og(n) is such that P(W n (rj, z) = 0) 
is the term of a convergent series. 



On logarithmic fluctuations. Our approach is inspired by Lawler, Bramson and Grif- 
feath's work jlOj . We develop three original ideas: (i) we propose a cluster growth model, the 
flashing processes, whose covering mechanism is simpler than internal DLA, (ii) we look at 
an intermediary scale, the scale of cells, since the deviations of the number of visits decrease 
with the cell-length, (hi) we build a coupling between internal DLA and flashing processes 
which allows to transport bounds from one model to the other. 

Let us describe how the idea of an intermediary scale is used in the context of flashing 
processes. Recall that we first partition Z d into a sequence of concentric shells. Each shell is 
partitioned into cells whose side length equals the width of the shell. Now, we observe that 
a site has good chances to lie inside the cluster if some cell, say C about this site, is crossed 
by many explorers. The notation W(n, C) refers to the number of explorers visiting C, when 
their initial configuration is n. We drop the index n appearing in W n (n, z) since there are no 
more constraint on not escaping the ball B(0,n). Now, the coupon-collector nature of the 
covering mechanism suggests that for some positive constant a^, 

W(rj,C) > a d \C\ x log(|C|) C C A(N) with a large probability. (1.13) 
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We neglect in these heuristics the log(|C|) term in (1.13). 

Note that in [TO], all the explorers start from the origin, whereas here, we only know 
that they cross C. For internal DLA, estimating the probability that C is not covered, when 
C is large and W(rj,C) > ctd|C| raises a difficulty which is absent when considering flashing 
processes. 

We now make our argument more precise. For a scale h and an integer K > 1, to be 
determined, assume that B(0, n — Kh) is covered by settled explorers. Partition the shell 
S = 1(0, n - {K - l)/i)\B(0,n - Kh) into about (n/hy- 1 cells, each of volume h d . It is 
also convenient to stop the explorers as they reach the boundary of M(0,n — Kh). Thus, 
with such a stopped process, explorers are either settled inside B(0, n — Kh) or unsettled 
but stopped on its boundary, that we denote by <9B(0, n — Kh). What we have called earlier 
the number of explorers crossing C is taken here to be the unsettled explorers stopped on 
C n <9B(0, n - Kh). 



Assuming (1.13) holds, it remains to show that the probability of the event {3C G S : 



W(rj,C) < ctd\C\} is small. We improve (1.11) by first using the independence between 



W(rj,C) and M(B(0,n - Kh),C), and then by replacing A R {fb,ri) by 1(0, n - Kh) in flLlOfr 
with R = n- Kh and n = |B(0,n)| 1 , 

W(r], C) + M(B(0, n-Kh),C)> M(rj, C). (1.14) 

Also, we define 

H(C) = E[M(rj, C)\ - £[M(B(0, n - Kh),C)]. 

Now, using that M(r/,z) and M(B(0, n — Kh),z)) are sums of independent Bernoulli, we 
show that ( 1.14[ ) implies a gaussian type lower tail 



P ( W( V ,C) < a d \C\) < exp |^- (MC) c;/( ^ |C|r I , (1-15) 

for a positive constant c, and where v{C) 

v(C) = var (M(rj, C)) - var (M(B(0, n - Kh),C)) . 

We then show that both jj,(C) and v{C) are of order K\C\. Then, P(W(r],C) < ctd\C\) is 
summable as soon as K\C\ > Alog(n). 



Recent works on internal DLA. After this paper was submitted, Jerison, Levine and 
Sheffield pi] improved our bound on the outer error in d = 2, by following a different approach. 
Then, by refining our approach, we obtained in [1] a bound of order ^\og(n) for both internal 
and external errors in dimension three or more. Finally, Jerison, Levine and Sheffield [H] by 
using their approach proved similar bounds in dimension three or more. 



Outline of the paper The rest of the paper is organized as follows. Section [2] introduces 
the main notation, and recall well known useful facts. In Section |3j we build the flashing 
process, give an alternative construction through exploration waves, and sketch the proof of 
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Lemma L2 In Section |4j we prove Propositions 1A\ and |1.5 using the construction in terms 
of exploration waves. In Section [5j we obtain a sharp estimate on the expected number of 
explorers crossing a given cell, and prove feature (ii) of the flashing times. Both proofs are 
based on classical potential theory estimates. Finally, in the Appendix, we give a proof of 
Lemma 



1.2, and recall a result of Sebastien Blachere. 



2 Notation and useful tools 
2.1 Notation 

We say that z, z' G Z d are nearest neighbors when \\z — z'\\ = 1, and we write z ~ z'. For 
any subset A C Z d , we define 

<9A = {z G Z d \A : 3z' eA,z'~z}, 

For any r < R we define the annulus 

A(r,R) = B(0,R) \B(0,r) and A(r, R) = A(r, R) n 1 d (2.1) 

A trajectory S is a discrete nearest-neighbor path on Z d . That is S : N — > Z d with S'(t) ~ 
S(t + 1) for all £. For a subset A in Z d , and a trajectory S, we define the hitting time of A as 

H(A; S) = min{t > : S(t) G A}. 

We often omit S in the notation when no confusion is possible. We use the shorthand 
notation 

B n = B(0,n), B n = B(0,ra), H R = H(B C R ), and H z = H({z}). 

For any a, b in K we write a A b = min{a,6}, and a V b = max{a, b}. Let T be a finite 
collection of trajectories on Z, d . For R > 0, z in Z d and A a subset of Z d , we call M(r, R, z) 
(resp. M(T,R,A)) the number of trajectories which exit B(0,i?) on z (resp. in A): 

M(T,R,z)=J2 1 {S(H R )=zh and M(r, J R,A) = ^M(r, J R,z). 

ser zeA 

When we deal with a collection of independent random trajectories, we rather specify its 
initial configuration rj G N z , so that M(r), R, z) is the number of random walks starting 
from 7] and hitting B(0, R) c on z. Two types of initial configurations are important here: (i) 
the configuration nl z * formed by n walkers starting on a given site z*, (ii) for A C Z d , the 
configuration 1a that we simply identify with A. For any configuration rj G N zd we write 



\V\ 

For any A C Z d , we define Green's function restricted to A, G\, as follows. For x,y G A, 
the expectation with respect to the law of the simple random walk started at x, is denoted 
with Eje, (the law is denoted ¥ x ) and 



G A (x,y) = E x 



1 {S(n)=y} 

0<n<H(A c ) 
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In dimension 3 or more, Green's function on the whole space is well defined and denoted G. 
That is, for any x, y G Z d 



G(x,y)=E x 

In dimension 2, the potential kernel plays the role of Green's function 



1 is(n)=y} 

,n>0 



a(x, y) = lim E a 



J2(l{S(l) = x}-l{S(l)=y}) 



1=0 



2.2 Some useful tools 

We recall here some well known facts. Some of them are proved for the reader's convenience. 
This section can be skipped at a first reading. 

In [TU], the authors emphasized the fact that the spherical limiting shape of internal 
DLA was intimately linked to strong isotropy properties of Green's function. This isotropy 
is expressed by the following asymptotics (Theorem 4.3.1 of [H]). In d > 3, there is a 
constant K g , such that for any z ^ 0, 



G(0,z) 



G, 



I d-2 



< 



with Cd 



(2.2) 



v d (d-2)' 

where Vd stands for the volume of the euclidean unit ball in M. d . The first order expansion 



(2.2) is proved in [TTJ for general symmetric walks with finite d + 3 moments and vanishing 



third moment. All the estimates we use are eventually based on (2.2) and we emphasize 



the fact that the estimate is uniform in ||z||. There is a similar expansion for the potential 
kernel. Theorem 4.4.4. of [11] establishes that for z 0, (with 7 the Euler constant) 



fn \ 2 i (\\ in 2 ^ + lQ g( 8 ) 
a(0, z) - - log ( z ) 

7T 7T 



< 



12 ' 



(2.3) 



We recall a rough but useful result about the exit site distribution from a sphere. This is 
Lemma 1.7.4 of [gj. 

Lemma 2.1 There are two positive constants C\,C2 such that for any z G dB(0, n), and 
n > 



Cl 



rr 



T < n(s(H n 



z) < 



C-2 



d-l ' 



(2.4) 



We now state an elementary lemma. 

Lemma 2.2 Each z* in 7L d \ {0} has a nearest-neighbor z (i.e. z* ~ z) such that 
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Proof. Without loss of generality we can assume that all the coordinates of z* are non- 
negative. Let us denote by b the maximum of these coordinates and note that 

\\z*\\ 2 < db 2 , and b > 1. (2.6) 

Denote by z the nearest-neighbor obtained from z* by decreasing by one unit a maximum 
coordinate. Using (2.6) 

||^|| 2 -||^|| 2 = 6 2 -(6-l) 2 = 26-l>6> t^. (2.7) 

yd 



Note that Q follows from 2||z*||(||z*|| - ||z||) > ||-2*|| 2 - \\z\\ 2 , and fl2J[ ). I 
We state now a handy estimate dealing with sums of independent Bernoulli variables. 

Lemma 2.3 Let {X n , Y n , n G N} be independent 0-1 Bernoulli variables. For integers n, m 
MS = Xl H \-X n and S' =Y x -\ h Y m . Define fort el 

/(t) = e*-l-*, and #(t) = (e* - l) 2 . 

If0<t< log(2), then 

IKgigS * - > s - s '' + £ • <-> 

Assume now that for k > 1, swp n E\Y„\ < ^=1. Ift<0, then 



K-l 



(2.9) 



Proof. Let X be a Bernoulli variable, and p = i?[X]. Using the inequality e x > 1 + x for 
x G K, we have 

E [exp (t(X - £[X]))] = pe^- p) + (1 - p)e~ tp = e~ pt (l + p(e* - 1)) < exp {f{t)E[X}) . 

(2.10) 

For a lower bound, we distinguish two cases. 

First, assume t > 0. We claim that exp(x — x 2 ) < 1 + x for < x < 1. Indeed, we use 
three obvious inequalities: e x > 1 + x for x G R, (i) for x < 1, 1 + x + x 2 > e x , and (ii) 
(1 + x 2 )(l + x) > 1 + x + x 2 . Thus, 

e x \l + x) > (1 + x 2 )(l + x) > 1 + x + x 2 > e x . 

This yields the claim. Now, set x = p(e* — 1), so that x < 1 when e* < 2. The last inequality 
in ( gig} yields 

£ [exp (t(X - E[X]))] > exp {-tp + p(e* - 1) - p 2 (e t - l) 2 ) = e-^W*)** 2 . (2.11) 
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Assume now that t < 0, and for k > 1, p < - — . We claim that for < x < — 

' ' r ft K 

K 



exp ( —x — —x ) < 1 — x. 



(2.12) 



Indeed, we have an additional inequality (iii) 1 — x + y > exp(— x) when x > 0. Note also 
that 



1 + -x 2 (1 - x) > 1 - x + 



x 



X < 



K — 1 



Thus 



e^ x ll -x)> ( 1 + ^x 2 ) (1 - x) > 1 - x + ^ > e" x 



2 " J x ~ - ~2 

k-i 



Now, set x = — p(e l — 1) > 0, so that x < We obtain 

E [exp - > exp (-tp + pie 1 - 1) - ^ 2 (e* - l) 2 ) = e 



/(*)p-f<?(i)p 2 



(2.13) 



The inequalities (2.8) and (2.9) follow (2.11) and (2.13) 



3 The flashing process 

In this section, we construct the flashing process, and state the crucial "Uniform hitting 
property" . We then present a useful equivalent construction in terms of exploration waves. 



Finally, we explain the coupling of Lemma 1.2 and postpone its proof to the Appendix. 



3.1 Construction of the process 

Partitioning the lattice. We are given a sequence H = {h n , n G N}. We partition the 
lattice into shells (Sj : j > 0). For an illustration, see Figure [X} For a given parameter 
ho > the first shell So is the ball B(0, ho). For j > 1, shell j is the annulus (see its definition 



(2.D) 



S; 



{rj-h^rj + hj), 



where {rj, j > 1} is defined inductively by r% — ho + hi, and for j > 1 

O+i ~ hj+i = + hj. 

In Section El we need that (o) V, is increasing, (i) j h-» hj/rj is decreasing, and (ii) hj = 

1 /3 

0(r ■ ). These properties are a straightforward consequence of our hypothesis hj < hj + i < 
(1 + jr)hj. Actually we will only need these properties and our hypothesis is no more than 
a sufficient condition. 



We also define 



S = {0} and Sj 



,rj), 3 > 1. 
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Flashing times. The key feature we expect from the flashing process is that its covering 
mechanism be simple. More precisely, our construction is guided by property (ii) of the 
Introduction which states that the flashing position, in a shell, is essentially uniform over 
the cell an explorer first hits upon entering the shell. Thus, we need to define together cells 
and flashing times to realize property (ii). It is important that all sites of a shell can be 
chosen as flashing sites with about the same frequency. In this respect, let us remark that 
a cell in shell Sj cannot be a ball of radius hj centered on E 3 -. Indeed, if this were the case, 
sites at a distance about hj would be in much fewer cells, than sites of Ej, and this would 
fail to make the covering of a shell uniform. We find it convenient to build a cell with a 
mixture of balls and annuli. A (random) flag Yj tells the explorers whether it flashes upon 
exiting either a sphere or the boundary of an annulus, whose distance from Ej is governed 
with a random radius Rj of appropriate density. Also, to allow for the possibility of flashing 
on its hitting position on Ej, we introduce an additional flag Xj. 

More precisely, consider {Xj,Yj,j > 0} a sequence of independent Bernoulli variables 
such that 

P(X j = l) = l-P(X j = 0) 



hf 



and, 



p<y j = i) = i-p<y j = o) = { \ 



Consider also a sequence of continuous independent variables {Rj,j > 0} each of which has 
density gj : [0, hj] — > IR + with 

„, dh d ~ l . , 

9j(h) = -j-5". (3-1) 

For j > 0, and Zj in Ej, let S be a random walk starting in Zj, an define a stopping time a 
as follows. If Rj = h for some h < hj then 

ifXj = l, 

H(M(z j: h A (rj + hj - \\zj ||)) c ) if Xj = and Yj = 1, 
H{A(rj - h, r-j + hf) if Xj = and Yj = 0. 

We set Hj = H(Ej), and we define the stopping times (aj : j > 0) as 

Oj = Hj + a(S o 9 Hj ), 

where 9 stands for the usual time-shift operator. For j > we note that, by construction, 
S(t) G Sj for all t such that Hj < t < o~j and we say that o~j is a flashing time when S{<jj) 
is contained in the intersection between Sj and the cone with base B(S(Hj), hj/2). We call 
such an intersection a cell centered at S(Hj), that we denote C(S(Hj)). In other words, for 
any z G Ej 

C(z) = Sj n{x eR d : 3A > 0, 3y G B(z, hj/2),x = \y} . (3.2) 



The uniform hitting property The main property of the hitting time a constructed 
above is the following proposition, which yields property (ii) of the flashing process to be 
defined soon. 
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Proposition 3.1 There are two positive constants oci < 0:2, such that, for ho large enough, 
j > ; Zj G Ej, and z* G C(zj). 



!<p,(5 W -.)<! 



(3.3) 



The proof of Proposition 3.1 is given in Section (5} 



The flashing process. Consider a family of N independent random walks (Si : 1 < i < N) 
with their stopping times (-Hjj, <J%j : j > 0). Let also Zij = Si(Hij) be the first hitting 
position of Si on Ej. 

We define the cluster inductively. Set A*(0) = 0. For i > 1, we define r* as the first 
flashing time associated with 5j when the explorer stands outside A* (2 — 1). In other words, 

t* = min {a itj : j > 0, S^K,-) G C(*y) D A*(z — l) c } , 

and 

A*(i)=A*(i-l)U{5 i (7f)}. 



3.2 Exploration Waves 

Rather than building A*(N) following the whole journey of one explorer after another, we 
can build A* (N) as an increasing union of clusters formed by stopping explorers on successive 
shells. Similar wave constructions are introduced in [10] and [9]. We use this alternative 



construction in the proof of Propositions 1.4 and 1.5 

We denote by G (Z d ) N the explorers positions after the k-th wave. We denote by A* k (N) 
the set of sites where settled explorers are after the k-th wave. Our inductive construction 
will be such that 

£ E fe & G U j<k Sj & G A* k (N). 

For k = we set £ (i) = 0, and ^(i) = 0, for 1 < i < N. Assume that for k > 0, A* k (i) is 
built for % = 0, . . . , N. We set A* k+1 (0) = A* k (N). For i in {1, • • • , N}, we set the following. 

• If £fc(z) ^ E fe , then 

6+1 (*) = e u i<fe 5 i> and ^fc+i(*) = ^fc+i(^ _ !)■ 

• If f fc (i) G E fc and Sifa,*) G C(z i)jfc ) fl ^(z - l) c , then 

6k+i(i) = $(<7 i)fc ) € £fc> and -^fc+iW = - !) u {^(^fc)} • 



• If G E fc and ^(a ijfc ) £ C(zj, fc ) fl *4£(i - l) c , then 

6+1 (*) = Si(H itk+1 ) G Efc+i, and A* k+1 (i) = A* k+l (i - 1). 
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In words, for each k > 1, during the k-th wave of exploration, the unsettled explorers move 
one after the other in the order of their labels until either settling in Sf--i, or reaching 
where they stop. We then define A*(N) by 

A*(N) = \jAl(N). 

k>l 

We explain now why this construction yields the same cluster as our previous definition. 
An explorer cannot settle inside a shell it has left, and thus cannot settle in any shell Sj 
with j < k if it reaches Now, since each wave of exploration is organized according to 
the label ordering, the fact that an explorer has to wait for the following explorers before 
proceeding its journey beyond does not interfere with the site where it eventually settles. 



3.3 Coupling internal DLA and flashing processes 



Proof of Lemma 1.2 For each positive integer N, we build a coupling between A(N) and 
A*(N). We first describe the main features of our coupling in words. Its precise definition 
is postponed to the Appendix. 

We launch N independent random walks, and build inductively the associated clusters 
A(l), A(2), . . . , A(N). In doing so, we use the increments of these random walks to define, 
step by step, N flashing trajectories S^, . . . , up to some times t\, . . . , i^. Let us describe 
informally step i + 1 of the induction. Assume that SI, . . . , S* are defined up to some 
times t% < ti, . . . , ^ < ti, and that each site of A(i) is covered by exactly one Sl(t k ) with 
1 < k < i. We can think of Sl(ti), . . . , S*(ti) as the positions of stopped flashing explorers, 
some of them stopped at one of their flashing times - say on blue sites - some of them 
not - say on red sites. Then, we add the i + 1-th explorer and flashing explorer. We set 
*S7+i(0) = <Si+i(0) = 0- We add new increments both to S i+ i and to the trajectory of one 
flashing explorer, say with label j in {1; . . . ; i + 1}, in such a way that the current position of 
the walker i + 1 and that of the flashing explorer j coincide. The label j is defined inductively 
as follows. Initially, j = % + 1. Assume that the walker i + 1 flashes on a red or blue site 
inside A{i). This site is occupied by exactly two stopped flashing explorers, j and j' (and all 
other red and blue sites of A{i) are occupied by exactly one flashing explorer). Since flashing 
explorers can settle at their flashing times, it makes sense, when j is flashing, to add the 
next increment to the trajectory of flashing explorer j 1 rather than j. We do so in two cases. 
First, when this happens on a red site. In this case, we turn blue that site since j is stopped 
at a flashing time. Second, when this happens on a blue site, say z, and j' > j. Note that in 
this case, both explorers flash on z but explorer j reaches z before explorer j' when launched 
in their label order. Our choice is such that the eventual cluster A* (N) has the correct law. 
In all other cases, we keep adding the increments of S^+i to the same flashing trajectory. It 
is important to note that the value of the increment does not depend on the index of the 
trajectory we choose to extend. Walker i + 1 eventually steps outside A{i). When the last 
walker steps outside A(N — 1) we have 

A{N) = {S* 1 {t 1 );---;S* N (t N )} with \A(N)\=N. (3.4) 

To define A*(N) we launch again, in their labels order, our flashing explorers from their 
current positions (possibly some or none of them since some or all of them can already have 
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reached their settling position). We then get 

A*(N) = {SKt*);--- ;S* n (t* n )} with \A*(N)\=N and r k > t k , for all k. (3.5) 



Proof of Corollary 1.3, Since a flashing explorer that visited some site beyond a given 
shell cannot settle in that shell, the one-to-one map 

i> N :S* k (t k )eA(N)^S* k (T*)eA*(N), k = l,...,N (3.6) 

satisfies, for all k and /, 

Sl(i k )^\JS m =► Sl(T* k )=^ N (Sl(i k ))?\JS m . (3.7) 



m<l 



m<l 



Thus, for all N > there is a coupling and a one-to-one map ipN between A(N) and 
A*(N) such that for all k > 1 

r/> N {A{N)nW rk+hk ) d A*{N)nW rk+hk . 



The inclusion has two important consequences. 

a. If A*(N) C M rk+hk , then A(N) C !,-»+/*. Indeed, any site in A(N) outside B n+hk 
produces, through ip N , a site in A*(iV) outside M rk+hk . 

b. If B rfc+/ife C A* (AT), then B rk+hk C A(AT). Indeed, those sites in A(N) that are mapped 
through ipN on A*(N) nM rk+ h k = M rk+ h k are necessarily contained in B^^. Since 
their number is |B rfc+?tfe | and ipN is one-to-one, they completely cover M rk+hk . 



4 Fluctuations 



In this section we prove Propositions 1.4 and 1.5 To do so we use the construction in 



terms of exploration waves of Section 3.2, Thus, we think of the growing cluster as evolving 
in discrete time, where time counts the number of exploration waves. The proofs in this 
section rely on potential theory estimates which we have gathered in Section |5j for the ease 
of reading. 



4.1 Tiles 



We recall that we have defined a cell of Sj in (3.2), as the intersection of a cone with Sj. We 
need also a smaller shape. We define, for any Zj in Sj, and for a small eo to be defined later, 



C{zj) = Sj n{x eR d : 3\>0,3ye B{ Zj , e hj),x = Xy} . 



(4.1) 
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As in Lemma 12 in [9], concerning locally finite coverings, we claim that, for ho large enough, 
there exist a positive constant Kp, and, for each j > 0, a subset Ej of Ej such that 



Wy G 5,- 



j^GSj-: ye CO) J < AT F and 5,-= (J Cfo) 



(4.2) 



For any Zj G Ej, we call tile centered at Zj, the intersections of C(zj) with Ej. We denote by 
T(zj) a tile centered at Zj, and by Tj the set of tiles associated with the shell Sj\ 

(4.3) 



(4.4) 



7~ = < T{zj) : zj e E, 



We choose eo to satisfy two properties. First, for any z e 5 3 -, there is Zj G Ej such that 

ze f| C(y). 

2/eT(%) 



This is ensured by the choice of a small enough eo- Indeed, let Zj G Ej be a site realizing the 
minimum of {\\z — y\\ : ?/ G Ej}. There is A > and u G B(zj, 1), such that 2; = Am. Now, 
there is Zj G Ej such that \\zj — Zj\\ < eohj, and for any y G T(%), we have ||y — Zj|| < 2e$hj. 
Thus, for eo small enough so that 1 + 2e$hj < hj/2 



g n~ Zj ) 



\u — y\\ < \\u — Z j 



\zj — y\\ < 1 + 2e hj < 



fl; 



j > 1, and any tile T G Tj 



which implies (|4.4|). Secondly, the size of a tile should be such that for some k > 1, for any 

At — 1 , 

(4.5) 



sup P z (S(H(Zj)) G T) < 



The inequality (4.5) follows from Lemma 5(b) of [10] (or Lemma 5.1 below) which for a 
constant Jd yields 

sup P 2 (S(H(Zj)) G T) < JrfJJL 

*eB(Q,r J — ftj-) 



The choice of eo is such that Jd\T\ < ^rhj 



4.2 Bounding inner fluctuations 

For n > 0, we take N = |B n |, we recall that *4|(iV) C ^ +1 (iV) for k G N, and *4*(iV) = 
U fc >!^(A^). We consider 

T* = min j fc > 1 : (J 5,- t A* k (N) L . (4.6) 

Note that A%(N) C Uj<fc5j, so that T* is the first time k when the At th wave does not cover 
all its allowed space. We recall that time is intended as counting the number of exploration 
waves. 
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For the flashing process if Uj<kSj (£ A* k (N), then for any k! > k, we have Uj^Sj <jL 
A* k ,(N), so that T* is also the shell label where the first hole of A*(N) appear. We have, for 
/ with r t < n, 



P(T* <l) = P (1(0, r t + h{) t A*{ N)) < P ( T * = k + l). 



(4.7) 



k<l 



In this section, we estimate from above the probability P(T* — k + 1) assuming r k < n. 
For k > 1 and A C we call Wk(A) the number of unsettled explorers that stand in A 



after the k-th wave, that is 



N 



W fc (A) = ^l A 



(4.8) 



i=i 



We now look at the crossings of tiles of Tk- On the one hand, we will use that if Wk(T) is 
large, then it is unlikely that a hole appears in the cell containing T during the k + l th -wave. 
We use for this purpose the fact that covering for the flashing process is similar to filling 
an album for a coupon-collector model. On the other hand, if r k is small it is unlikely that 
Wfc(7~) be small. We now make precise what we intend by small and large. For any positive 
constant £, we write 



P (T* = k + 1) =P (T* = k + 1, VT e T k , W k {T) > o 

+ P(T* = A ; + l,3TeT fe ,W fe (T) <0 



<P (t* = k + i|vr g T fe , W k (T) >t)+p(3Te Tk, W k (T) < e) . 



(4.9) 



A coupon-collector estimate. The first term in the right hand side of (4.9) is bounded 



using a simple coupon-collector argument. Indeed, the event {T* = k + 1} implies that there 



is an uncovered site in S k , say z, when explorers stopped in E fc are released. By (4.4), there 
is z k G Sfc, such that z is a possible settling position of all explorers stopped in T(z k ). Now, 
knowing that {Wk(T(zk}) > £}, Proposition 3.1 tells us that the probability of not covering 
this site is less than (1 — ai/hf) to the power £. In other words, 

p (t* = k + i|vr g r fe , w fc cn > e) < |5 fc | ^i - j^J < \s k \ ex P • 



Henceforth, we set 



£ = 74/i d log(ra) with fa = sup{h k : r fc < n} 



fc:r^,<n 



(4.10) 



and A large enough so that 

P(T* = k + l\\/TeT k ,W k (T) >£) < |B n |exp(-aiAlogn) < ^. (4.11) 
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Estimating {W k (T) < £}. For any T G 7fc, we consider the counting variable Lk(T) 
M(B(0,r fe — h k ),r k ,T), and define 



M k (T) = W k (T) + M(A* k ,r k ,T), so that M k (T) = w M(iVl {0} , r k , T). 



(4.12) 



The idea of defining M k and (for the internal DLA process), and bounding W k by M k — L k , 
is introduced in [10J. Our main observation is that L k (J~) is independent of W k (T), and 

W k (T) + L k (T)>M k (T). 

As a consequence, for any positive constants t and £ (and with the notation X = X — E[X]), 

P (W k (T) < <e^ x E [exp (-tW k (T))} 

_ t z E[exp(-t(W k (T) + L k (T)))] 



E[ex V (-tL k (T))] 
< exp (-t (E [M k (T) - L k (T)} - 0) x 



E[exp(-t (M fc (T)))] 
£ [exp (-tl k (T))] 



Using Lemma 2.3 with condition (4.5), we obtain 



log P (W k (T) <£)<-t(E [M fc (T) - L*(r)] - + /H)£ [M fe (T) - L k (T)} 

+ ^9(~t) Py(S(H(i: k ))eT). 



rk-h h ) 



We now proceed in two steps. We show in Step 1 that for some constant k' 



E [M fc (T) - L fc (T)] > «V - (r k - h k 



h 



d-l 



„d-l ■ 



Since {h k /r k , k > 0} is non increasing, it follows that there is a constant K\ > such that, 
for all a > 0, and fc a := sup{j GN: rj < n — a/tlogn}, where ft, is defined in (4.10), we 
have ^ 

inf £ [M fc (T) - ^fc(T)] > «/(n d - (n - h) d )^— > Kl ah d logn. 

k<k a n» 1 



Now, if we choose £ as in (4.10), with a = 2A/k\ and A;* = fc Q , i.e., 



A;* := sup jj G N : r 3 - < n ftlog(n) j 



then, we get, for all k < k*, 

E[M k (T)-L k (T)]>2^. 
We show in Step 2, that for a constant C depending on the dimension only 

K (S(H(Z k )) eT) < C E [M fc (T) - L k {T)\ . 



(4.13) 



(4.14) 
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Suppose for a moment that Step 1 and Step 2 hold. Since, for some c > 0, max(/(- 1), g(— t)) < 
ct 2 when t < 1, there is d > such that for k < k* 



logP (W k (T) < Ah d log{n)) < inf f-t + c(l + ^)t 2 ) £ [M fc (T) - L fc (T)] 

< - c'E [M fc (r) - L fc (r)] < -2d Ah d log(n). 
Now, using (4.9), (4.11) and ( 4.15[ ) for A large enough, we have 



(4.15) 



P( T * = k)< 



k<k* 



n- 



Borel-Cantelli's lemma yields then the inner control of Proposition 1.4 



Step 1. We invoke Corollary 5.4, with n = r k , and A n = h k (the hypotheses h k = 0(r^ ) 
and /ifc large enough hold here, as seen in the first paragraph of Section 3.1). We have for 
some positive constants k', K, and for n large enough 

E[M k {T)-L k {T)) =E[M((\M n \ - |B rfc _ h J)l ,r fc ,r)] 

+ £?[M(|B rfc _ h Jl , r k , T)\ - E[M(B rk _ hk ,r k , T)} 



> 



^ hk \)¥ (S(H k )eT)-Kh 



d-i 



>2n\n d - (r k - h k 



,d\ l0 k 



d-l 1Xn k 



d-1 



>«V-(^-M d )TT. 

r k 



for r k < n and Hq large enough. 



(4.16) 



Step 2. By Lemma 5.1 below, there is a constant kq such that, for y G B(0,rfc — h k ), and 

z et k 

F y (S(H(Z k ))eT(z))< ,lV;|T(: - )l 



\z - y 



\d-l 



Therefore 



P 2 y (S(H(Z k ))eT(z))< E 



J 



2(d-l) 



yeB(0,r fc -h fc ) j:h k <j<2r k y:j<\z-y\<j+l 

For a constant C^, we bound |{y : < \z — y\ < k + 1} < Cdk d ~ l . Thus, 
J2 K(S(H(Z k ))eT(z))< J2 CdK % iZ) l2 



y&(0,r k -h k ) 



j:h k <j<2r k 



<C'\T{z)\ 2 f l d=2 log(n) + l d>2 -^ 



(4.17) 



(4.18) 



Since \T(z)\ is of order hfr 1 , (1414] holds 
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4.3 Bounding outer fluctuations 

This section follows closely [9j. The features of the flashing process allow for some simpli- 
fication. We keep the notation of the previous subsection. There, we proved that for some 
integer k*, which depends on n, 

P(T* > k*) = l-e(n), with ^e(n)<+oo. 

n>l 

The integer k* is the largest such that r k * < n — 2Ah log(n)//ti, for a large constant A and 
with h defined in (4.10). As a consequence, the following conditional law can be seen as a 
slight modification of P. 

P*(.) = P(.\T* > k*) (4.19) 

We begin by proving that under P* the probability to find some k with n < r k < 2n and 
some tile T in T k with W k (T) larger than or equal to £' = 2A'h d logra for a large enough A' 
decreases faster than any given power of n. First, note that on {T* > k*} 

W k (T) + Li(T)<M k (T), with L* k = M(M(0,r k *-h k *),r k ,T). (4.20) 

Our key observation is that the pair (W k (T), l{T*>k*}) is independent of L* k . Thus, for any 
t > 0, 

P(W k (T) > e, T* > k*) <e-*'E [e tw ^ I {T *>** } ] 

te E[exp(t(W k (T) + Lp) l {T * >k » } ] 
E [e tL l] 

E [e tL l] 

= exp (-t (C - E [M fc (T) - LI))) x 1 J . 

£/ le fcj 

By Lemma 2.3, we have 
logP (W h (T) > e, T*>k*)<-t (£' - E [M fc (T) - L*]) + /(f) x J5 [M fc (T) - L*] 

j/eB(0,r fc .-/n..) 

The steps are now similar to the previous proof. We first estimate E [M k (T) — L* k \. By 
Corollary |5.4[ for some positive constant K' and for n large enough, 

E [M k (T) - L* k (T)] < K' (n d - (r>) d ) ^ + 0(h d k ^) 

"Vr ( 4 - 22 ) 
< K'dn d -\n - r fe .)-fcr + 0(^ 1 ). 
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Note that since r k < 2n, we have r 



d-i 



o[n 



i-i, 



n — r k *)) so that 0(h k 



is small compared 



to the first term in (4.22). Since k H- h k /r k is decreasing, we have for some constant K 

E [M k (T) - L* k (T)) < Kh d \ogn. 



Secondly, we estimate the sum of P 2 (S(H(E k )) G T) which appears on (4.21 ). We use (4.17) 
again to obtain as in (4.18) and for a constant C 



E 

,rk*-hk*) 



F 2 y (s(H(i: k ))er)<chl id ^ 



l d=2 log(n) + l d>2 



(r k - 0> - h k *)) c 



Note that r k — (r k * — h k *) > h k , and, since k H- h k /r k is decreasing, we have, for n large 
enough, h k < h k *{r k /r k *) < h x (2n)/(n/2). Thus, for a constant C, 



£ {S{H{E k )) eT)<C( l d=2 h 2 log(n) + I d>2 /i d ) < C/i d log( 



We choose A' = K to obtain 

log P WT) > T* >k*)<- (Kt - Kf(t) - Cg(t)) ^log(n). 

Since we have P(T* > k*) > 1/2 for A large enough and K can be taken as large as we 
want, we have that P*{W k {T) > £') decreases faster than any given power of n. 

Now, let F k denote the event that no tile T in E^ contains more than £' = 2A'h d logn 
unsettled explorers after the fc-th exploration wave. We define, with the notation of Section 
[3j Qk = c(£o, . . . ,£*.), and note that F k and {T* > k*} are (^-measurable. 

For any tile T G %, let G Sfc be such that T = T(z k ) and denote by C = C(z k ). We 
are entitled, by Proposition 3.1 to use a coupon-collector estimate on the number of settled 
explorers during the k + 1-th exploration wave. On F k PI {T* > k*}, and for some positive 
constant K\ 



E 



\A* k+1 nC\ 



> \c 



> \c 




-Qti 



w fc (r) ] 
^ J 



> KxW h {T) inf 

x<2A'logn 



X 



We now write for some positive constant K 2 



inf 

x<2A' logn 



> 



inf 



j g— a\x/2A' logn 



> 



2^4'logn x<2A'io g n x/2y4'logn 

I 1 _ e -ais ^ ^ 

x ~ logn 



. inf 

2,4' logn x<i 
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We conclude that on F k n {T* > k*}, 



E 



A* k+1 nC 



>K X K, 



W k (T) 
logn 



(4.23) 



Recall now that property (4.2) implies that Kp\A* k+1 PI S k \ > Yl z e£ l^fc+i ^ C( z k)\- Thus 



summing over ^ G with C = C(z k ) and T = T(z k ) in (4.23), we obtain on F k f]{T* > k*} 



E 



\A* k+1 n 5, 



> if 



W fc (5 fc 



where X 



Also, since Wfc(«Sfc) < 



logn 

n)\, we have, for n large enough, 



^ LlF fc n{T*>fc*}|A+i H^fclJ > K — — n f{^ k ). 



Since P(T* > fc*) > 1/2, 

£*[|^ +1 nS fc |] > ^*P"2n^). 
In other words, noting that \A* k+1 PI <Sfc| = W k (S k ) — W k+ i(S k+ \) 



E* [W k+1 (S k+1 )} < 1 



K 
logn 



E* [W k {S k )} + 2n d P{F c k ). 



(4.24) 



(4.25) 



By iterating ( 4.25[ ), and using our previous estimate on P*(W k (T) > £'), we obtain that for a 
large enough e, P*[^ n+<E iog 2 n(<^ n +eiog 2 n)]5 * s summable, when l n is the lowest index for which 
r i n > n - Also, the probability (under P !) of seeing at least one explorer reaching the shell 



S, 



i n +elog 2 n 



is summable. Using Borel-Cantelli lemma, this yields the proof of Proposition 1.4 



4.4 Lower bound for the deviations 



4.4.1 Proof of Proposition 1.5: the outer deviation 



We denote by K n the largest index such that Sk„ C B(0, n), and by E n the event that all 
explorers stopped on T*k„, at time K n , settle afterward in one of the shells {Sj : K n < j < 
K n + blog(n)} for some positive constant b, and note that E n = {A*(N) C Uj<K n +biog(n) Sj}. 

We want to find b such that Yln>i E(E n ) < oo. Using that the flashing times of the 
different explorers are independent, we have 

P(E n ) <E[P (all explorers, stopped in Y, Kn , flash in U j<Kn +bio g (n) Sj\Q Kn )] 



<E[( sup P(an explorer, started on z, flashes in Uj<x„+fciog(n) Sj\Gx n )) 



]■ 



Also, there are at least \Sk„\ explorers stopped on T,x n , and there is a positive eo such that 
the probability of crossing a given shell without flashing is larger than eo- Thus, 

|*-*iC I 

P{E n ) <E [( 1 — inf P(an explorer started on z is unsettled at time K n + 61og(n)) ) 



< 1 



felog(n) 
-0 



When choosing 6 small enough, we reach 5^ n >i E(E n ) < oo. 
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4.4.2 Proof of Proposition [1.5^ the inner deviation 



We recall that K n is the largest index such that Sk k C B(0, n). The rough idea here is that 
when we stop explorers on S^ n _i there are necessarily tiles (of Ek„-i) containing of the 
order of h d ^ 1 _ l sites and which receive h d K _ x explorers. The number of explorers on these 
tiles is not enough to cover the associated cells with the coupon collector mechanism. We 
now make rigorous such an argument for shells with index of order K n — log(K n ). 

To simplify the notation, let us first define three positive constants ci, c% and C3 such that 
for any k with n/2 < r k < n, we have 



\<Sk\ < cih k ri 



d-1 



K 



sup z£S , 



'z,6h k ) fl Ei 



|Efc| > C2n d l h k , and inf \C(z)\ > c^h 

z£T, k 



Using «2 given in Proposition |3.1[ we define 

1 



8a 2 



\og(h Kn ), and A n 



C2C3 



(4.26) 
(4.27) 



Now we assume ho large enough to have A n a strictly positive integer. 

We wish now to consider a peel of A n shells before <9B(0,n). Let I n be the index of 
the inner shell in this peel, that is r In+ A n < n < r In+ A n +i- Since T* < I n + 1 implies 
that Uj<i n <Sj (£ A*(N), it is enough to show that P(T* > I n + 1) decays faster than any 
polynomial in n. 

Note that the monotonicity of k h-» h k /r k and rj n > n/2, imply that 2hi n > h k for 
I n < k < I n + A n , and n large enough. Also, on the event {T* > I n + 1}, we have 
,n n — hi n ) = A* In (N) after the I n -th wave. Thus, 



W In (E In ) = |B(0,n)| - |B(0,r /n - < 2c 1 A n n d - i x fc Jn . 

A key feature of the flashing process is that explorers stopped, at time /„, outside M(z, 3h In )n 
Ej n cannot settle in C(z). In other words, knowing Q k , the covering of a family of cells 
{C(zj),j = 1, . . . ,Af} are independent events if \\Zi — Zj\\ > 6hj n for i 7^ j. Now, there is an 
integer Af and sites {zj, j = 1, . . . , Af} with 



Vz ^ j 



> 6h In , and ^ m Zj ,6h In ) nE/J > |E/J. 



We then get using (4.26) 



Let 



in — 2 



d-1 



(4.28) 



r = { 3 e[i,Af}:W In (B(z J ,3h In )nj: In )<c 3 a n h d I j, and r c = [i,A/]\r. 

On {T* >I n + 1}, 

2c 1 A n n d ~ 1 x h In > W In (E In ) > W In (B( Zj , 3h In ) n E 7 J > |r c | x (c 3 a n /il) . 
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Thus, using the definition of A n in (4.27), and the bound (4.28) on Af, we obtain 



< 2dA n h In n d 1 < c x c 2 c z a n M = M_ 
~ c 3 a n h d In ~ 2c 3 ci a n c 2 2 ' 

In other words, we have that |T| > Af/2. Now, as already noticed, knowing Qj n , for any 
subset / C [1,A/], the events {C(zj) C A* In+1 (N), j 6 /} are independent. By conditioning 
on Qj n , we obtain for ho large enough 



P{{T* > I n + 1}) =E 



-E 



\I\>N/2 



lr=/ xP Vie J.Cfo) c X +1 (N)\g It 



/C[1,AH, UI>W/2 is/ 



< sup P(A* In+1 (N)DC( Zj ), W In (B( Zj ,3h In )ni: In )<c 3 h d In a r 



N/2 



(4.29) 



Considering the probability appearing on the right hand side of (4.29), we can think of a 



coupon-collector problem, where an album of size \C(z)\ has to be filled when we collect no 
more than c^hj a n coupons. Using inequality (4.31) of Lemma 



4.1 



below, we show that 



P ({T* > I n + 1}) < exp (-^a^l^-V- 1 



This concludes the proof. 

The result about filling an album, that we just mentioned, is based on the following 



simple coupon-collector lemma (together with Proposition 3.1 ), which we did not find in the 
vast literature on such problems. 

Lemma 4.1 Consider an album of L items for which are bought independent random cou- 
pons, each of them covering one (or possibly none) of the possible L items. IfYi is the item 
associated with the i-th coupons, we assume that for positive constants ai,a 2 , such that for 
any j = 1, . . . , L, 

rt/i A'r, 

(4.30) 



Let tl be the number of coupons needed to complete the album. Then, for any < A < 
4^1og(L) ; we have 



P{t l < AL) < exp 



a 



2j^2 e ~2a2A 



L < exp 



a\A 2 



(4.31) 



Proof. We denote by Oi the time needed to collect the i-th distinct item after having 
collected i — 1 distinct items. The sequence {a\, a 2 , ■ ■ ■ , ctl} is not independent, but if 
34 = cr({Yi, . . . , Y k }), and r(k) = G\ H h cr k , then for i = 1, . . . , L 



ai(L — i + 



>p{(T i >k\y T(i . 1) )> i 



a 2 (L — i + 



(4.32) 
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Indeed, calling E{% — 1) the set of the first i — 1 collected items, 

P (cri> k y T (i-i)j =P f{K r (i-i)+i, • • • ,Y T (i-i) + k} C £(i - 1) y T (i-i) 



P [Y e£(i-t 



X-(i-l) 



1) 



Using (4.30) we deduce (4.32) from (4.33). Formula (4.32) gives that 



ai(L — z + lj a 2 (-^ — 2 + 1; 



as well as 



E 



y T (i-i) 



< 2 



af(L-z + l) 2 ' 



Now, we look for 5 < \/Z such that 



(4.33) 



(4.34) 



BVL 



Note that 



E i a L-i] > 2AL - 

i=y/L 

e^]>-e4t^ Ai °^) 

/ — 1. QJ5 2 + 1 tt9 



(4.35) 



Thus, condition (4.35) holds for £> > exp(2a2A), but recall that B < yL also, and this gives 
a bound on A. Finally, note that 



max I E 



cr L -i 



y r(L — i— 1) 



and set 



X: 



E 




yr(L-i-l) 









L,...,BVLl < 



< 1. 



Vl 



a i 



For x < 1, note that e x < 1 + x + x 2 to obtain for < A < 1, by successive conditioning 

(BVL \ / By/L 

i=^/L ) \i=vT 

< e -A Ql AVZ "Q j^ 1 + A 2 supE [Xfl^^^n < 1-30 

i=sfh 



< exp I -XaiAVL + A 2 ^ sup E X 2 y r{L 



-i-i) 
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Finally, we have, using (4.34) 



^2 sn P E x f y T {L-i-i) < ^2 a i su p — 

i=\/L i=\/~L 



a 



L-i 



The results follows as we optimize on A < 1 in the upper bound in (4.36). 



5 Potential theory estimates 



We collect in this section three technical results. In Corollary 5^, we estimate the differ- 
ence between the expected number of independent random walks exiting a ball B(0,n) at a 
distinguished site, whether the random walks are initially on the origin or are spread over 
a sphere B(0,r n ) with r n < n. Corollary 5.4 is used to bound the mean number of explor- 



ers exiting some large ball from a given site, and its proof relies on a discrete mean value 



property Theorem |5.2[ which in turns relies on Blachere's Proposition |B.1| written in the 
Appendix. Then, Lemma |5.1| improves an estimate of Lawler, Bramson and Griffeath in 



[TU] dealing with the exit site distribution from a sphere when the initial position is not the 
origin. Indeed, Lemma 5(b) of [TU], states that when d > 2, there is a positive constant 



such that for any r > 0, z G B(0, r), and z* G 



W z (S(H r 



< 



r), we have 
J, 



\d-l 



(5.1) 



Thus, when \\z*\\ — \\z\\ is small, ( 5.1[ ) is useless. Since we need bounds on the sum of squares 
of F z (S(H r ) = z*) over z G B(0,r — h) of order log(r) in d = 2, and of order l/h d ~ 2 when 
d > 2, we establish the following. 

Lemma 5.1 There is a positive constant kq such that, for all r > 0, if z G B r , z* G dM r , 
then 

F z (S(H r ) = z*) < „„ (5.2) 



\z — z 



Finally, we prove the uniform hitting property, Proposition 3A, for the boundary of a cell. 
Though the property is natural, the non-spherical nature of a cell, makes its proof tedious. 



5.1 A discrete mean value theorem 

The following result has interest on its own. 

Theorem 5.2 There are positive constants Kq and K a such that for any sequence {A n ,n G 
N} with K < A n < n 1 / 3 , for any z G B n with n — \\z\\ < I, we have, setting r n = n — A n; 

| |B r J x G n (0,z) - G n(y,z) I < K a . (5.3) 

y& rn 
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Remark 5.3 Note that a related (but distinct) property was also at the heart of [10J. 
Namely, for e > 0, and n large enough, if z G B n , and n — ||z|| > en, 



l n \xG n (0,z) > ^G n (y,z). 



(5.4) 



We start with proving the following useful corollary of Theorem 5.2 



Corollary 5.4 In the setting of Theorem \5.£\ and for any A C dM n 

E[M(\M rn \l ,n,A)}- E[M(M rn ,n,A)}\ < K a \A\. 

Proof. Note that (5.5) holds if for any z* G <9B n 

||B,J x F (S(H n ) = z*)- ^y(S(H n ) = z*)\< K a . 



(5.5) 



(5.6) 



By a classical decomposition (Lemma 6.3.6 of [II]), we have for a finite subset B C Z d , 
y £ B, and 2* G 

F v (S(H(dB) = z*) = ± G ^*)- ( 5 - 7 ) 



For i? = B(0, n), we replace in (5.6) the value of ¥ y (S(H(dB)) = z*) by the right hand side 
in (5.7), and are left with proving that for any z6l„ with n — \\z\\ < 1, we have (5.3). I 

Proof of Theorem |5.2[ When d > 3, we express G n (0, z) in term of Green's function 
(Proposition 4.6.2(a) of [IT]). 

G n (0, z) = G(0, z) - E,[G(0, S(H n ))}. 



Now, using Green's function asymptotics (2.2), there is a constant K\ (independent on n) 
such that 



v d G n (0,z)-2°^\ 

n a-i 



< K }, where a(z) = E z [\\S(H n 



(5. 



In d = 2, G n is expressed in terms of the potential kernel (Proposition 4.6.2(b) of [TT] ) 

G n {0,z) = -a(0,z)+E z [a(0,S(H n ))]. 



Using (2.3), we have 



7rG n (0, z) = 2a(z)/n + 0{l/n 2 ). 



d-i 



Now, r d = n d — dA n n' 
0(n 1/3 ), and < n - \\z\\ < 1 



0(A„n ), so that using (5.8), and the hypothesis A n 



G„(0, a! ) = (rS + O(r*- 1 ))(2^g + O(l) 



(7/ - dA^ 1 + C(A^ d - 2 ) + O^- 1 )) + O(^)) 

2a(z){n - dA n ) + O(l) 



(5.9) 
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Since {HS'nll 2 — n, n £ N} is a martingale (with the natural filtration), 

E z [\\S(H n )\\ 2 } - \\z\\ 2 } = E z [H n ] = G n (y,z). 

Using n — ||2;|| < 1, this yields for a constant Ki 



} J G n (y,z) - 2a(z)n 



<K, 



(5.10) 



We now invoke Proposition B.l of the Appendix. There is Kb such that for z £ B n with 
n — \\z\\ < 1 



^ G n (y, z) - 2a (2;)«iA r 

2/e«4(r ra ,n) 



< .K";,, where 0:0(2) = E 2 



\\S{H n )\\-\\z\\ H n <H{B rn 

(5.11) 



From (5.10) and (5.11), we obtain 



} J G n (y, z) - 2na(z) + 2a (z)dA n 
y m rn 



<K, + K b 



(5.12) 



Now, from (5.9) and (5.12) we obtain for a constant K 2 , 



„\G n (0,z)- G n(y,z)\ -2(a (z)-a(z))dA r 



Now, from 

\a (z)-a(z)\ <¥ z (H(B rn ) < H n ) x (a (z) +E Z [\\S(H n )\\ 
and the gambler's ruin estimate, for K Q > and z £ A(n — 1, n) 

F z (H(B rn )<H n ) < 

we deduce that 



H n >H(B r S\) 



A ' 



A n \a (z) - a{z)\ < 2A n F z (H{B rn ) < H n ) < 2K . 
The desired result follows. 

5.2 Proof of Lemma [5JJ. 

We follow the proof of Lemma 5(b) of pU]. Set D := \\z — z*\\. Let O' be a closest point to 
(1 + §)z* in M(z*, f ). We define := B(0', § ), 5^ := B(0', f ), and we note that 

\\z-z*\\ < \\z-0'\\ 
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and, for all x in dB' 2 , the triangle inequality \\z — z*\\ < \\z — x\\ + \\x — z*\\ implies that 

„ D 

min \\z — x\\ > — . 

x£dJ3' 2 3 

Now, define 

t := inf {t > : S(t) G {z} U 5 r c } , and r := inf {i > : S(t) e B[U dB' 2 } . 

By a last exit decomposition together with the strong Markov property 

F z (S(H r ) = z*) =G r {z,z)F z *{S{r) = z) 

<G r (z,z)FAS{r') E B' 2 ) max FJS(t) = z) 

x£dB' 2 

=F Z * (S(r') e B') max G r (x, z) ( 5 - 13 ) 

x£dB' 2 

<F z *(S(t) G B') max G r+D (x,z). 

x£dB' 2 

A gambler's ruin estimate yields, for some positive constant c 

FAS{r')eB' 2 )<^. 



The desired result follows from (|5.13|) and the previous bound after we show that for a 

D 
3 



constant c, such that for all x satisfying ||x — z\\ > D 



G r+D (x,z)<^- 2 (5.14) 

On the set V := B(z, f ), the map y i— > G(x,y) is harmonic. By Harnack's inequality, we 
have 

G r+D &z) < ^dJ2 G r+n&y) = (5-15) 

where c is a positive constant and Y is the number of visits of V before time H r+ D- By 
taking the supremum over the entering site of V in (5.15) 



G r+D (x,z) < — supEjF]. 

L> y&V 



It remains to show that sup y&v F, y [Y] < JD 2 , for some positive constant J. This is identical 
to (2.10) of [10J, and we omit this last step. 



5.3 Proof of Proposition 3.1 



For j > 0, consider Zj in Sj. We show that for positive constants a±, «2, and for all z* in 
C(zj) we have (3.3). The random walk has initial condition S(0) = Zj. 



First, when z* 



"31 



S{<Tj 



if and only if X~ = 1. This happens with probability 



1/hj, and gives the result in this case. 
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Assume z* G C(zj)\{zj}. We recall that the unbiased Bernoulli variable Yj decides 
whether the explorer can flash upon exiting either a sphere or an annulus. More precisely, 
we draw Rj with density gj given in (3.1), and if Yj = 1 (resp. Yj = 0) the walk flashes 



upon exiting the ball of center Zj and radius Rj A (r,- + hj — \\zj\\) (resp. 



Rj^j + Rj)) 



provided S(tjj) G C(zj). 

Step 1: Flashing when exiting a sphere (Yj = 1). We first prove the upper bound 



when Yj = 1 and Xj = 0. It is obvious that 



z* G 9E(zj, \\z* — Zj\\), but z* £dB(zj,\\z* - Zj\\ - 1). 
Thus, Rj G]|| Z II ' 1 1 

111, and there is a constant C such that 



P(Xj = 0, ^ = l,Rj e]\\z* - Zj\\ - 1, ||** - Zj\\\) < C 



\Z — Zi 



\d-l 



hj 



(5.16) 



On the other hand, by (2.4) of Section 2.2 
P (S{<jj) = 



X, = 0, Yj = 1, Rj e]\\z* - Zj\\ - 1, \\z* - Zj 



< 



c 2 



\Z* — Zn 



\d-l 



(5.17) 



The upper bound in the case {Xj =0, Yj — 1} follows from (5.16) and (5.17). 

We turn now to the lower bound when Yj = 1 and Xj = 0. Since we want a lower bound, 
we consider the event that the walk flashes on z* when exiting a sphere only in the case 
where | — Tj\ < hj/2. Note that by Lemma 2.2, z* has a nearest neighbor, say z, which 
satisfies 

II _ .11 < II *_ .1' 1 

1 1 Z Z j 1 1 1 1 Z Z j i 



This means that if h G V :- 



Zj\\ - \\z* - Zj\\[, then z* G dM(zj, h). Thus, 



P z . {S{<jj) = z*) > P{Xj = 0, Yj = 1, Rj eV)x inf P 2j (S(H(dB( Zj , h))) 



> n ^mf v F Z3 (S(H(dM( Zj ,h))) = z* 

ch d ' x ci 
> , x 



hj h d ~ l 



(using (2.4)). 



The lower bound in the case Yj = l,Xj = 0, and |||**|| — Tj\ < hj/2 is obtained. 

Step 2: Flashing when exiting an annulus (Yj = 0). The upper bound for this case is 
close to the case Yl- = 1. It is obvious that 



z* G dA{rj — |||**|| — Tj\, Tj+\ 



\z — r 



|), but ** 9A(rj— |||2;*||-rj|+l,rj + |||2;*||-r^|-l). 



Thus, necessarily i?j G]|||**|| — Tj\ — 1, |||**|| — Tj\], and 

P(15 = 0,i2, e] I || - Tj\ - 1,|||**|| - r,|]) < d 
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For h > 0, define T>h = A(rj — h, rj + h). It is enough to prove that for some constant c, and 
for any h such that z* G <92\, (and h G]|||z*|| — Tj\ — 1, |||z*|| — Tj|]) 

P*, (^M) = **) < j^n- (5-18) 

Note the following fact. If ||z*|| > \\zj\\, and the walk exits T>h at z* , then the walk exits 
j,Tj + h) at z*, whereas if \\z*\\ < \\zj\\, and the walk exits T>h at z*, then the walk enters 
j,Tj — h) at z*. In both cases, Lemma 5(b) of [10] yields then (5.18). (Actually, Lemma 
5(b) of [TO] is formulated to cover only the case \\z*\\ > \\zj\\, but its proof covers both cases). 



We turn now to the lower bound. By Lemma 2.2 z* has a nearest neighbor, say z, 

. 1 . 

|||z|| -rj\ < \ \\z || -r,j\ - — . (5.19) 

This means that if h e V := [\\\z*\\ - r 5 \ - l/(4y/d), \ \\z*\\ - rj\[, then z* e &D h . 

We only need to consider the case \\\z*\\ — rA > hj/2. It is enough to prove, for h E V, 
z* G C(zj) PI dV h , and for some constant c (that depends on d) 

^ (S(H(V c h )) = z*) > (5.20) 

Let y* be the closest site of <9B(0, rj) to the segment [0, z*], and let x* be in M rf given by 

/ z* 
X = ^ + 2 ] M- 



Note that if z ~ z and z satisfies (5.19), then ||x* — z\\ < inf 2e gx) h \\x* — and we define 

R* = -(\\x* -z\\ + inf 
Define 

^ = A(r,-|, ri + |), mds et r = i (l ',ir)n«. 

Thus, if \\z*\\ > \\zj\\, then T is the boundary of the lower hemisphere of the ball M(x*,R*). 
We need also the time r + = inf{n > 1 : S(n) G T> c h U {^j}}. By a last exit decomposition, 
and the strong Markov property, we have 

P,. (S(H(&D h )) = z*) = G Vh (z v z,)W z * (S(t + ) = Zj ) 

> G Vh {z vZj )V z , (H(T) < r+)mmF x (S(r) = z 3 ) (5 _ 21) 



>P,. (H(T)<r + ) mm G Vh (x, 



.7/ 



Since z* G C(zj), we have \\y* — Zj\\ < hj/2, so that y* and Zj can be connected by 10 
overlapping balls of radius hj /10 in such a way that, applying Harnack's inequality 10 times 
(see Theorem 6.3.9 in [11]) to the harmonic map y i— > Gx> h {x, y), we can estimate from below 
the last factor in (5.21). For any x G T 

G Vh (x,Zj) > c^G Vh (x,y*). 
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We use again Harnack's inequality on the harmonic functions x i-> Gx> k (x,y*), to obtain 

mmG Vh (x,y*) > c H G Vh (x', y*), 



where x' G B(x*,i?*/2) and — G [| — 1, |]. The purpose of choosing x' is to have 
y* G B(x',h/A), and B(V, /i/2) C V h so that G Vh {x',y*) > G M{x , >h/2) {x',y*). 

When dimension is 2, the classical expansion of GM( x ',h/2)( x ' \ •) ( see Proposition 6.3.5 of 
[TT] ) gives with a constant K 2 



Gn(x',h/2)(x',y*) > -lo. 



7T 



/i/2 



r^^.og2-l^. (,22) 



When /i is large enough, G B{x ^ h/2 )(x' , y*) > \og(2)/n. 

When dimension is larger than 2, by using (2.2), there is a constant Kj, such that, when 
h is large enough 



G Vh (x',y*) =G(x',y*) - E x > [G{S{H{V c h )),y*)] 
> 



h d-2 \ hj ~ hd ~ r 



(5.23) 



As a consequence of (5.23), we just need to prove that the first factor in (5.21) is of order 
1/h at least. We realize the event {H(T) < r + } in two moves: the walk first hits the sphere 
M(x*, R*/2), and then exits from the cap dM(x*, R*)(lV h . 



F z * (H(T) < r+) >^P, (H(B(x*,R*/2)) < H(M c (x* , R*))) 



x 



inf F y (S(H(B c (x*,R*))) G V h ) . 

<x*,R*/2) y V V ' J 



(5.24) 



The first factor in the right hand side of (5.24) is of order 1/R*, that is of order 1/h. To 
deal with the second factor, we invoke Harnack's inequality to have for y G dB(x*, R*/2), 
and for x" the closest point of 7L d to x*, 



F y (S(H(B c (x*,R*))) G Vf) > c H F x „ (S(H(M c (x*,R*))) G V, 



(5.25) 



We invoke now (2.4) to obtain for some constant K3 



F x » (S(H(M c (x*,R*))) eVh) > Cl 



\dM(x*,R*)nV h \ 



\dB(x*,R* 



(5.26) 



We gather (5.24), (5.25) and (5.26) to obtain the desired lower bound. 



APPENDIX 
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A Coupling and proof of Lemma \1.2 



We give a precise definition of our coupling. To avoid heavy notation we write the coupling 
algorithm pseudo-code. 

First of all, we draw N independent sequences of independent Bernoulli and continuous 
random variables ((Xk t i,Ykj, Rk,i '■ I > 0) : 1 < k < N) as in Section [3] In addition, we call 
{Uk '■ k > 1) the sequence of the increments of a generic independent simple random walk 
on Z d . From these two sources of randomness we extract our explorer and flashing explorer 
trajectories with their associated clusters. The flashing times will be adapted to the flashing 
explorer trajectories as in Section [3j 

Our pseudo-code is made of two loops of size N that make precise the previous description. 
With the first loop we build our N random walk trajectories ((Si(t) : < t < ri) : 1 < i < N) 
with their associated clusters A(l), . . . , A(N). Step by step, within this first loop, we also 
define pieces of the flashing explorers trajectories S^, . . . , S^. With the second loop we 
complete the trajectories of the flashing explorers to build the associated cluster A*(N). 
During the algorithm, G N stands for the time up to which the trajectory of flashing 
explorer k has been defined (k G {1; ■ ■ • ; N}). We use the same t for the time governing 
the evolution of each simple random walk S{. The index j is updated before adding each 
random walk increment to the partial sum of Si and S*. The updating procedure uses the 



index j' described in Section 3^, and we denote by A = U the increment. Each encountered 
U stands for the first unused random variable in the sequence (Uk '■ k > 1). 

The main advantage of the pseudo-code formalism is that it allows, through the as- 
signment operator expressions of the kind j t— max(j, j') or tj 4— t j + 1 rather than 
j(9+l) = max(j(0), j'(9)) and tjf$ + x\(9+l) = tj^(9) + l with 9 a discrete parameter ordering 
the sequence of our elementary moves. It makes also implicit identities like tk(9 + 1) = 
for any quantity tf. that does not need to be updated. Our following pseudo-code can be re- 
written in a classical inductive way with 9 running through {(z, t) G {!;••• ; iV} x N : t < r^} 
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according to lexicographic order. Marks ^ and refer to remarks (a) and (b) below. 

A(0) G- 0; 
For % = 1 to iV 

j <- i; 

t G- 0; tj G- 0; 
Si(t)<-0; S*(^)g-0; 
While Si(t) G A(i — 1) 

If tj is a flashing time for explorer j( a \ then 
unique index^ k e {!;■■■ ;i}\{j} 



(Note that Si(t) = SUtj) ) 



3 <~ 



such that Sl(t k ) = SUtj) = $(*); 



If tji is not a flashing time for explorer j'^ a \ then j G- j'; 
Otherwise j G- max(j,/); 
A G- [/; + 1) G- Si(t) + A; + 1) G- + A; 

(so tftoi S^t + 1) = Sjfo + 1); 

t<-t + l; tj<-tj + l; 
<- A(i - 1) U {$(*)}; 
k i G- i + 1; 

A*(0) G-0; 
For fc = 1 to N 

While tk is not a flashing time for explorer or S%(tk) G ^4fc_i 
3S(t fc + l)<-SjE(t fc ) + tf; 

A*(fc)G-A*(fc'-l)U{S*(* fe )}; 
fc G- fc + 1; 



Remarks: 



(a) Recall that for / = j,j' or k, Sf is defined up to time t\ as well as its associated flashing 
times. 

(b) One checks by induction on i that just after the instruction u A(i) G- A(i — 1) U {Si(t)}" 
we have 

A(i) = {^(0;---;5?(*i)} and L4(i)|=i (A.l) 
To do so, one checks by induction on t < r, h that 

A{i - 1) = {Slit,)- ■■■■ S^it^)- S* +1 (t j+1 ); S*(t t )} and \A(i - 1)| = i - 1. 

(A.2) 

Since we always have S*(tj) = Si(t) this proves by induction that f is well defined. 

The key observation is that for each increment U, the index of the explorer that follows 
this increment depends on the whole previous construction, but the value of U does not 
depend on it. As a consequence, we build independent random walks S±, . . . , Sn coupled 
with independent flashing random walks S^, . . . , S^. Then, one simply checks by induction 
on % and k that 



A(i) = {5i(Ti);...;5 i (r i )} and A*(k) = {S^t*); • • • ; ^(r fc *)} 



(A.3) 
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for all 1 < i, k < N. 
Finally, define 



, tiv) an d ( r i\ • • • , t n) the values of (ti, . . . , tjv) at the end of the 



first and last cycle respectively. Since t\ 



tjsr can only increase during our loops, we have 



r* > t k for all k. Then (3.4) and (3.5) follow from (A.l) and (A.3). 



B Time spent in an annulus (By S.Blachere) 

This section is devoted to an asymptotic expansion of the expected time spent in an annulus 
A(r n ,n) for r n < n, when the random walk is started at some point z within the annulus, 
and before it exits the outer shell. 

Proposition B.l There are positive constants Kq, Kb, such that for any sequence {A n , n G 
N} with Kq < A n < n 1 ^ 3 , for any z G A(r n , n), we have setting r n = n — A n; 



with 



G n (z,y) - {2dA a {z) - d( 

yGA(r n ,n) 

a (z) = E z \\\S(H n )\\- 



"-Ikll) 2 ) 



<K b ((n-\\z\\)Vl) , 



(b.i; 



H(B c (0,n))<H(B(0,r n )) 



Proof. Our strategy is to decompose a path into successive strands lying entirely in the 
annulus. The first strand is special since the starting point is any z G A(r n ,n). The other 
strands, if any, start all on dM(0,r n ). We estimate the time spent inside the annulus for 
each strand. Let us remark that we make use of three facts: (i) precise asymptotics for 
Green's function, (ii) ((7(0, S(n)), n G N) is a martingale, and (iii) ( || /S'(rz) || 2 — n, n G N) is 
a martingale. 

Choose z G A(r n ,n). We define the following stopping times (Dj, Ui, i > 0), correspond- 
ing to the i th downward and upward crossings of the sphere of radius r n . Let 9(ri) act on 
trajectories by time-translation of ra-units. Let r = H(B Tn ) A H n , D = Uq = 0, and 

Dx = rl HiBrn)<Hn + ool Hn<H{Brn) . 

If D\ < oo, then U\ = H Tn o 9(Dx) + D\, whereas if D\ = oo, then we set U\ = oo. We now 
proceed by induction, and assume D i: Ui are defined. If Di = oo, then D i+1 = oo, whereas 
if Di < oo, (and necessarily Ui < oo) then 

A+i = Ui + (rl T=H{Brn) + ool T=Hn ) o 6(Ui), and U i+1 = D i+1 + H rn a 9(D i+1 ). 

With this notation, we can write 



G n (z, y) = E z [r] +^E 2 [ro 9(Ui)l Di 

y&A{r n ,n) i=l 

= E z [t]+F z pi < oo) x 1(2;), 



<OOj 



(B.2) 
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where 



oo i— 1 

l(z) = ^E 2 [r o 0(1^)1 A < oo] J] (l - P z (l> i+1 = oo 

4=1 " j=l 



£>j < oo 



(B-3) 



Now, we compute each term of the right hand side of (B.2). 

We have divided the proof in three steps. 

Step 1: First, we show that there is a positive constant K, (independent of z and n) such 
that when z G A(r n ,n), then 



P, (A < OO) 



a (z) 



<^((n-P||)Vl) 



(B.4) 



Note that when z G B(0,n), and n — ||z|| < 1, (B.4) yields 

E 



P* {Di < oo) 



l|S(r)||-|M 



D 1 = oo 



A, 



< 



A 2 "' 

n 



(B.5) 



Secondly, we show that for z G A(r n ,n), and z > 1 



P, (a+i = oo 



Di < oo - 



E-. 



(\\S(U t )\\-\\S(D 



i+l)\\) - L L> 1 o0(C/ i )<oo 



Dj < oo 



< 



K 



A 2 ' 

n 



(B.6) 

Our starting point is the classical Gambler's ruin estimate, which in dimension 2 reads with 
the potential kernel instead of Green's function, 



P 2 (Dx < oo) 





G(0,z 


)-E 2 


G(0,5(r)) 


D\ = oo 




E z 


G(0,S(r)) 


Di < oo 




G(0,S(r)) 


Di = oo 



(B.7) 



We now expand Green's function (resp. the potential kernel) using asymptotics (2.2) (resp. 

define a 

(l|5(r)|| : 



(2.3)). For this purpose, it is convenient to define a random variable 

1 



X z) 



2 ) 



Note that for any z G A(r n ,n), X(z)/\\z\\ is small. Indeed, 

X(z) (\\S(r)\\-\\z\\)(\\S(r)\\ + \\z\ 



(B.8) 



Since A n = n — r n = Oi^n 1 ^ 3 ), we have for n large enough 



X W < 2(n + 1 )A„ < 8A._ ^ ffflV < ^ x 1 



(n - A r 



n 



z£A(r n ,n) V \\ Z 



n n z 



(B.9) 
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More precisely, X(z) is of order 2(||S'(r)|| — \\z\\). Indeed, A^ < n, and (B.8) yields 

'(\\S(r)\\-\\z\\r 



X(z) = 2(\\S(r)\\-\\z\\) + 



X(z)-2(\\S(t)\\-\\z\\) < 



1 

(B.10) 

When dimension d > 2, we set 77(d) = d ~Y L . In order to use Green's function asymptotics 



(2.2), we express S(t) in terms of X(z) as follows 

1 



\\S(r)\\ d - 2 \\z\\ d - 2 



1 ( X(z)^- V(d) 
1 + 



\z\ 



(B.ll) 



We have a constant _?Q such that 



1 



2 _ V ^ 



< 



71- 



(B.12) 



For d > 2 and any 2^0, (2.2), (B.ll) and (B.12) yields 



G(0,S(r)) - G(0, ,) - ,(d)C, I + " W + 1 X(Z)2 



INI** -1 ' 2 ||z|| 
In dimension 2, the potential kernel asymptotic yields for K 2 > 



< 



a(0, 5(r)) - a(0, z) 



1 /X(z) lX{z) 



71 \ \\z\\ 2 \\z\ 



< 



K, 



n- 



(B.13) 



(B.14) 



l/vr, and ffiT= 1/2. 



In view of (B.14), we assume henceforth that (B.13) holds, but in d = 2 we think of n(d)Cd 
2 



Using (B.7) and (B.13), we obtain 



(A < 00) 



where 





'X(z) 


D\ = 00 


-C(z 




E z 


'X(z) 


Di = 00 


-E z 




Dx < 00 


+ C(z)-C(z) + 0(±) 



(B.15) 



X 2 (.) 



.Di = 00 



and g(z ) = g^+l E . 



X 2 (^) 



Di < 00 



Using ( |B.9[ ), we have some rough estimates on C and C_. For any z € ,A(r n ,n), 

1 



C(z) = O 



A 



A, 



and fiW = 0('^)=0(£ 



(B.16) 
(B.17) 



Using (B.10), we have better estimates for C and C_. 



C(z) = ^-",f ll)2 + O ( (n - |NI)Vr ) , C(z) = d { - " ;: " - ' 



n 



O 



\z \\-r n VI 



n 



(B.18) 
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The rough estimates (B.17) together with (B.10) allow us to derive from (B.15) an estimate 
for F z (Di < oo), for any z G A(r n ,n). 

F z (Di < oo 



E z 


~\\S(t)\\ -\\z\\ 


D\ = oo 






E z 


)\S(r)\\-\\z\\ 


D\ = oo 


-E 2 


IISMii - 1| 2 || 


Di < oo 





(B.19) 



This yields p4| ) since a (z) < 1 + (n - ||*||) V 1 < 2(n - ||z||) V 1. 
Case where 2 G dB(0,r n ). 
On {Di = 00}, we have 

X(z) = 2(\\S(t)\\-\\z\\) + 0(- 1 



On {D 1 < 00}, we have 
This implies 



X(z) = 2(\\S(r)\\-\\z\\) + 0(-). 

n 



(B.20) 



C(z) = d^ + 0( — ), and C(z) = 0(i). 
\\z\\ n n 



(B.21) 



Thus, 

P 2 (D 1 = 00) 



2E 2 


NI-IWII 


D x < 00 


+ <2(*) + o(±) 


E z 


'x(z) 


D\ = 00 






Di < 00 


+ ^)-C( 2 ) + 0(i) 



E, 



N|-||5(r)|| 



D x < 00 



(B.22) 



A n + 0(1) 



E, 



N|-||5(r)|| 



Di < 00 



A, 



AH 



In order to obtain (B.6), we write (B.22) on {Di < 00}, and z = S(Ui) as follows. There is 
a constant K such that on the event {Di < 00}, 

Es W) [(115(^)11 -||5(r)||)l Dl0 ^ i)<00 ] 



Es^) [l£> j+ i=oo] - 



A n X P s(l/i) (Di < 00) 



< 



K 

aI 



(B.23) 



Note that ( |B.22| ) implies that P s(l/i ) (£>i < 00) = 1 + 0(1/A„), so that QB.23P reads as we 
integrate over {Z)j < 00} with respect to E 2 

E 2 [l A<00 (ll^H-ll^r)!!)!^^] 



P 2 (A+i = 00, Di < 00) 



A T 



< 



K¥ z {Di < 00) 



AH 



(B.24) 
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We obtain ( |B.6| ) as we divide both sides of ( |B.24fl by P 2 (A < oo). 
Step 2: We show now that for any z G w4(r n , n) we have 

|E* [r] - (dA n a Q (z) - d(n - \\z\\) 2 ) \ < K ((n - ||z||) V 1) . 



When z £ B n and n — ||z|| < 1,(B.25) reads 



|E 2 [r] - (dA n a (z) - d(n - \\z\\) 2 ) < K. 
When z G A(r n , n), and i > 1, we show that 



(B.25) 



(B.26) 



E, 



.D, < oo 



E, 



(115^)11-115(^)11) 



_D, < OO 



dA 2 . 



< 



a! 



(B.27) 

Using that {(^(n)!! 2 — n, n G N} is a martingale, and the optional sampling theorem (see 
Lemma 3 of [TU] ) 

E 2 [r] =E 2 [||5(r)|| 2 ] - \\z\\ 2 = ||*|| x E 2 [X{z)\ 



\z\\ x E 



X(z) L>i = oo ¥ z (Di = oo) + E 2 X(z) D x < oo P 2 (.Di < oo)) . 



Thus, using (B.15), simple algebra yields 



E 2 [t] = \\z\\ x ((C(z) - C{z))W z {D x < oo) + C{z)) + 0(1). 



(B.28) 



By recalling (B.18) and (B.4) 



E 2 [r] =d ( ((||,|| - r n f - (n - ||,||) 2 + 0(A n )) ('^ + 0( - (n 11:11 1 V ' 



An v A^ 
+ d{n- \\z\\) 2 + 0{{n- \\z\\) VI) 
=d(2||*|| - n - r„> (^) + d(n - \\z\\) 2 + 0((n - ||^||) V 1) 
--dA n a (z) - d(n - \\z\\) 2 + 0((n- \\z\\) V 1) 



This yields (B.25) 



Assume now that z G dB(0,r n ). From (B.28), we have 

E 2 [t] = \\z\\ x ((C(z) - C(z))W,{D x = oo) + C{z)) + 0(1). 



We use (B.6) and (B.21) to obtain 



E z [r] =\\z\\ (dA 2 n + 0(A n )) 



- \\S(r)\\ 



D x < oo 



=dA n E 2 



N|-||5(r)|| 



L>i < oo 



A„ 



+ 0(1). 



A 2 ' 

n 



0(1) 



(B.29) 
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Now, write (B.29) as follows. There is a constant K such that for any z G <9B(0,r n ) 



T 



E z [||z|| - ||5(r)||l Dl<00 ] 



dAl 



A n P a (£)i < oo) 



< 



K 

ah 



(B.30) 



Note that by ( |R22| ) A n P z (Di < oo) = A n + 0(1) and |||z|| - ||5(r)||l Dl<00 | < 1, thus 

E z [r] E x [\\z\\-\\S(t)\\1 Di<o0 ] 



dAl 



A n 



< 



K 
A~2 



(B.31) 



We replace z by in (B.31) under the event {Di < oo} to obtain 



E s{Ui) [r} E S{Ui) [(\\SmW - \\S(D l o6(U i ))\\)l Iht>e( p i)<00 ] 



dAl 



< 



K 
AH 



(B.32) 



We multiply both sides of (B.32) by 1d i; <oo, take the expectation on both side of (B.32), and 

divide by P Z (A < oo) to obtain ( |B.27| ). 

Step 3: For i > 1, we show the following bounds 



2 > 7» > 7= , where 7» = E 2 



(115^)11-115(1) 



-D, < oo 



(B.33) 



The upper bound is obvious. For the lower bound, first we restrict to {Di < oo}, so 
that Ui < oo. By Lemma 2.2, S(Uj) has a nearest neighbor x, within B(0,r ra ) such that 
l|S(A)|| - \\x\\ > 1/(2 VS), and flB.33p is immediate. 
Step 4: We show ( |B.l ) using (B.2). For p such that 1 < p < oo, let 



P i— 1 

o v = E * [t ° < oo] J] (l - 

i=l 3=1 



_i = oo 



.D j < oo 



(B.34) 



Now, (B.3) reads \(z) = linip^oo a p (this is the limit of an increasing sequence). We establish 
in this step that, for some constant K, any integer n 



lim 

p— >oo 



1 _ _fp_ 

dAl 



< 



K 
A? 



(B.35) 



Once we prove (B.35), we have all the bounds to estimate the right hand side of (B.2). 
Indeed, using (B.25), (B.4) and (B.35), we have 

E z [t] + F z (D 1 < oo) x \(z) =dA n a (z) - d(n - \\z\\f + O {{n - \\z\\) V 1) 



< z ) , ^ A" -INI) vi 



A r 



AH 



x 



(dAl + 0(A n )) 



2dA n a (z) -2d{n- \\z\\) 2 + O ((n - \\z\\) V 1) . 



(B.36) 
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In order now to prove (B.35), we introduce first some shorthand notation. For p and j 
positive integers 



1 - ocj = F z = oo Dj < ooj , and # 



E, 



T o d{U 3 ) 



Dj < oo 



With this notation (B.6) and (B.27) read as follows. 

K 



ot ; 



1L 



< — , and 



ft 



2K 



Let us rewrite (B.34) as 



< — , so that \aj-Pj\ < ^ 2 



p-i 



ai = 1 — /3i and a p = a p -i — P p J^[(l — ocj) for p > 1. 



(B.37) 
(B.38) 

(B.39) 



In order to establish (|B.35|), we show by induction that 

v 



with for p > 1 



3=1 
P-I 



T-r . , 2if 

e p = e p _i + —]_]_(!- a,) and e x 



j=i 



A? 



(B.40) 



(B.41) 



Note that it is easy to estimate e p from (B.41). By (B.38) and for Kq large enough there is 
a constant Ks such that 



2K 



1 + exp 

k=l 

2K ( ' 

n V fc=i 



E 

i=i 

E 11 

3=1 



2K 

2A n / / * AH^ 5 



2Kks 



Now, by (B.38) (B.40) holds for p = 1, and we assume it holds for p — 1. Then 

p— i p— i 

(i - - - < a? < a - /?p) n( i - «*) + 

i=i j=i 



(B.42) 



Then by (B.38), we have (B.40) with e p satisfying (B.41). 



Now (B.35) follows as we notice that Step 3 implies, together with (B.38) and for K t 



large enough, that 



lim ij(l- ai )=0. 



3=1 
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